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ABSTRACT 

A planet orbiting in a disk of planetesimals can experience an instability 
in which it migrates to smaller orbital radii. Resonant interactions between 
the planet and planetesimals remove angular momentum from the planetesimals, 
increasing their eccentricities. Subsequently, the planetesimals either collide with 
or are ejected by the planet, reducing the semimajor axis of the planet. If the 
surface density of planetesimals exceeds a critical value, corresponding to ~ 0.03 
solar masses of gas inside the orbit of Jupiter, the planet will migrate inward a 
large distance. This instability may explain the presence of Jupiter-mass objects 
in small orbits around nearby stars. 

In the standard theory of solar system formation, solid material orbiting in a gaseous disk 
accumulates to form small rocky or icy bodies called planetesimals ([i]). Protoplanets then 
form by accretion of planetesimals. If a protoplanet accretes roughly ten Earth masses 
(10M e ), it can then capture a gas envelope from the protoplanetary disk and become a gas 
giant like Jupiter. We now know of a number of Jupiter-mass objects orbiting solar-type 
stars well inside the radius where rocky material can condense (Table 1; ref. 2-8). For 
example, a planet orbits r Bootis at a distance of 0.0462 astronomical units (AU), where 
the equilibrium temperature of 1550K is higher than the condensation temperature of the 
most refractory minerals. Although some of these detections are controversial (^), we feel it 
is likely that most or all are real. It is difficult to understand how to form such planets in 
place. 

Although it is difficult to form planets at such small radii, once in place they can survive 
(\TT\). Thus it is natural to ask whether giant planets can form at orbital radii of a few AU 
and then migrate inward. One proposed migration mechanism involves the generation by 
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the planet of density waves in the gaseous protoplanetary disk, which cause the planet to 
spiral inward (|7Zj , |7I|) . The movement of the planet might be halted by short-range tidal 
or magnetic effects from the central star (17Z|) ; however, short-range stopping mechanisms 
cannot easily explain the objects in Table 1 with semimajor axes a p ^ 0.2 AU. 

Another migration scenario involves interactions between two or more Jupiter-mass 
planets fl7I| , [T3[ ). If two such planets form with a small enough separation their orbits are 
unstable, resulting in either a collision or an ejection. The typical outcome is a system with 
a massive planet in an eccentric orbit at 1-2 AU, similar to 70 Vir or HD 114762. In a 
tiny fraction of cases the periapsis of one of the planets will come close enough to the star 
that tidal evolution circularizes the orbit, resulting in a massive planet on a circular orbit 
inside 0.1 AU. As pointed out in ( |73D this low yield is difficult to reconcile with the observed 
frequency of such systems, roughly a few percent of nearby stars. A second difficulty is the 
low eccentricity and relatively large a p of the planet orbiting p Corona Borealis, which would 
require an additional damping mechanism. 

In this paper we discuss a migration mechanism in which the planet exchanges energy 
and angular momentum with the residual planetesimal disk through resonant gravitational 
interactions, gravitational scattering, and physical collisions. Long-distance migration is 
triggered (the "migration instability") if the surface density E of planetesimals exceeds a 
critical value E c derived below. The migration halts when E drops below E c , or when a 
significant fraction of scattered planetesimals strike the star. Since planetesimals cannot 
survive inside a few stellar radii (R*), the migration is certain to halt. If the mass of the 
planet is less than about three times the mass of Jupiter (Mj), the eccentricity of the planet's 
orbit is reduced by the migration; if not, the eccentricity may increase. 

The orbit of a planetesimal is specified by its semimajor axis a, eccentricity e, and 
inclination i; the corresponding elements for the planet are a p , e p , and i p . In our analytic 
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work, but not in out numerical work, we assume that the mutual inclination between the 
orbits of planetesimals and planet is <^ 0.1 radians. It will also be convenient to use the 
specific energy E = — GM*/(2a) and angular momentum L = [GM*a(l — e 2 )] 1 / 2 , where M* 
is the mass of the central star; we shall not distinguish between the total angular momentum 
L and its normal component Lcosi. We assume that the ratio \x = Mp/M* of planet mass 
to stellar mass is much less than unity = M Q /Mj ~ 1047 for Jupiter), as is the ratio 
of planetesimal mass to planet mass. 

Planet formation remains poorly understood, but we flank this obstacle by focusing on 
processes that occur after formation is nearly complete. We shall assume for simplicity that 
only a single massive planet forms. The newly formed planet is probably surrounded by 
an annular gap in the planetesimal disk, whose radial extent is ~ 2.5 times the Hill radius, 
h = ^a p (d). However, a planet's gravitational reach can exceed its Hill radius grasp by a 
larger factor, as we now explain. 

Consider a planetesimal with orbital period T = 2w(a 3 /GM^) 1 / 2 and e < 1, perturbed 
by an exterior planet with period T p > T, a p = a + Aa and e p ~ e. If the planetesimal is 
in a mean-motion resonance, so that (k + q)T ~ kT p , where k and q are positive integers, 
the planet can force chaotic perturbations of the planetesimal orbit. The semimajor axis a 
is nearly constant, since the planetesimal is in resonance, but the eccentricity e undergoes a 
random walk, gradually diffusing to larger values. Eventually the planetesimal crosses the 
planet's orbit, after a time fll^) 



where k G [0, q] and q^2. This time is given in units in which the gravitational constant G, 
M* + M p , and a p are equal to unity, so T p = 2tt. The minimum e at which the planetesimal 
can enter the Hill sphere of the planet is e m ~ (a p — h)/a — l (|73|). For a Jupiter-mass planet 
at a p = 5.2 AU, and a planetesimal in the 7:4 resonance with initial eccentricity eo = 0.05, 
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we find T R « 0.5 x 10 5 yr (0). 

The orbital phase of a planet-crossing planetesimal at times of conjunction with the 
planet is effectively random because of resonance overlap. Thus close encounters between 
the two objects are common. The first close encounter generally removes the planetesimal 
from the resonance, but leaves it in a planet-crossing orbit. A planet-crossing planetesimal 
undergoes a random walk in both e and a (or in L and E), but the Jacobi constant J = E — L 
is roughly constant flip. Scaling arguments suggest that the time to random walk out of 
the system should generally be between T p j \x and T p //i 2 , corresponding to a range 10 4 -10 7 
yr for Jupiter. This range is consistent with the median lifetime of Jupiter-family comets, 
5 x 10 5 yr (|73|), and with the lifetimes found in our numerical simulations. For example, if 
we start a planetesimal in the 7:4 resonance it evolves at roughly constant E or a while e 
random walks from 0.05 to ~ 0.3 over ~ 3.5 x 10 5 yr (Fig. 1). At this point it suffers a close 
encounter with Jupiter. It then follows curves of constant J until ejected after about 4 x 10 5 



Planetesimals may also suffer fates other than ejection: 

1. Collision with the star: If the planetesimal random walks to sufficiently small L (large 
e), it will be absorbed by the star. This channel can become important once the planet 
has migrated to a p <^ 5-10-R*. Planetesimals on resonant orbits with a <^ (a p — h)/2 
also tend to collide with the star, because these become highly eccentric before they 
become planet-crossing. 

2. Collision with the planet: A planetesimal may collide with the planet before being 
ejected. The probability is 1 — exp(— P), where 



yr. 




- 6- 



where p is the mean density of the planet. This result shows that about 20% of orbit- 
crossing planetesimals will strike Jupiter but most will be ejected (|7^). 

3. Long-term capture into mean-motion resonances: Planetesimals on planet-crossing 
orbits can be captured temporarily into resonances; this is the inverse of the process 
described by equation ([!]), which now gives the typical residence time in the resonance. 
The median residence time is generally short compared to the age of the solar system. 
However, a few planetesimals can be trapped for very long times near stable islands 
(§/3). Such trapping has been seen in numerical integrations of planet-crossing orbits 



Planetesimals start with negative energy and are ejected with positive energy; thus the 
ejection process must remove orbital energy from the planet, which moves closer to the star 
( \22 |) . This phenomenon is well-known in the context of our solar system: Fernandez and 
Ip fll5[ ) suggested that the ejection of planetesimals caused the orbit of Jupiter to shrink 
by 0.1-0.2 AU. This shrinkage can naturally explain the depletion of the outer asteroid belt 
( [7^1 , If the surface density of planetesimals is above the critical value S c , this process is 
unstable. To calculate S c , we let the semimajor axis of the planet shrink by Aa p , reducing 
its energy by 



AE p ~ 



GM*M p Aa p 



(3) 



Planetesimals newly captured into chaotic resonances will be removed from the system, either 
by ejection or by consumption by the planet or the star. The mass of planetesimals affected 
is 



AM ~ 2irp c aa p Yi(oia T 



(4) 



where a < 1 is a measure of the average a/a p of the affected planetesimals (\2dj) , and p c ~ 1 is 
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the capture probability. In disposing of a mass AM of planetesimals the planet loses energy 



where / is the fraction of the original planetesimal energy that is taken from the planet's 
orbit; we expect that / ~ 1 for a p ^> R* and becomes negative as a p — > R # , where the planet 
perturbs most planetesimals into the star. 

The migration process is unstable if |Ai? a | > |A^ P |, which requires that the density 
exceed the critical density S c : 



in other words if the mass in the planetesimal disk interior to the planet is of order M p , the 
planet can migrate nearly to the surface of the star. 

The migration halts when either (i) the local surface density of planetesimals falls below 
the critical value, or (ii) a significant fraction of the planetesimals plunge into the star. The 
local surface density is sure to fall to zero near the star. Solid bodies cannot condense at 
radii <^ 7i2*, and existing planetesimals whose orbits might evolve to smaller radii cannot 
long survive at distances ^ 2R Sf . The minimum semimajor axis achievable by the migration 
instability is thus a few to ten stellar radii or ~ 0.03-0.1 AU. 

We have simulated the evolution of a Jupiter-mass body in a planetesimal disk using 
the Opik approximation fllffi); see Fig. 2. The simulations assume that the surface density 
in planetesimals varies as S(r) = £©(1 AU/r) L5 , where £©, the surface density at 1 AU, is 
a free parameter. The total mass in the planetesimal disk within radius r is then 



and assuming solar metal abundance (Z = 0.02) the disk mass in gas interior to Jupiter's 
orbit is O.16M (£ /1O 3 g cm" 2 ). 




(5) 




(6) 



M D ~ 1.4 x 1O _3 M (£ /1O 3 g cm~ 2 )(r/ AU) 



0.5 



(7) 
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The energy transfer from planet to planetesimal is approximated as a succession of close 
encounters until the planetesimal is ejected, strikes the planet or strikes the star. Summing 
over many encounters, we derive the average efficiency of energy transfer from the planet 
to planetesimals in a given resonance. These efficiency factors are then used to calculate 
the evolution of the planet in a disk of a given initial mass profile. For the case where 
the planet interacts with the 2:1, 3:2, 2:3 and 1:2 resonances, as well as a broad resonance 
zone in the immediate vicinity of the planet, the onset of the migration instability is at 
S « 200 g cm" 2 (Fig. 2). Equation (|), evaluated with a p = 5.2 AU, a = 0.7, p c = 1, 
and / = 1 gives S c « 70 g cm -2 , corresponding to S « 160 g cm -2 , in good agreement 
with our numerical estimate. A much higher density S Q ~ 8000 g cm -2 is required for the 
instability to persist until the planet reaches 0.03 AU. 

In some cases a gas disk with the surface density required for the migration instability 
may be nearly gravitationally unstable. The criterion for local gravitational instability of a 
thin disk is Q = cfl/ (7rGS) < 1 where c is the sound speed and fl is the angular speed (Jj2% . 
Replacing £ by E c from equation @ we find 

Q = 12 (ttu) " 2 (sjtk) " 2 (m) S) " 2 (3 ■ (8) 

where T is the disk temperature. 

Ejecting a planetesimal removes energy AE and angular momentum AL from the planet. 
If \AE/AL\ > 1, then e p decreases, while e p increases if \AE/AL\ < 1. For a planet with 
M p <; 3Mj, |A£J/AL| > 1 because the planetesimal must approach within a few Hill radii 
of the planet to be ejected, and the Hill sphere does not extend into the region where 
\AE/AL\ < 1 when E m 0. In this case the ejection of small planetesimals tends to reduce 
e p and a p . However the planet is immersed in a bath of planetesimals, many of which have 
substantial masses. Interactions between the planet and these objects will tend to produce 
equipartition between the energy in radial motion of the planet and the planetesimals. If 
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most of the planetesimal mass is in objects of mass near m , the expected equilibrium 
eccentricity of the planet is e eq ~ (mo/Mp) 1 / 2 . This ranges from 0.06-0.006 for m between 
an Earth mass and a lunar mass. 

Seven objects in Table 1 have a p < 1 AU, and probably did not form in their present 
orbits. We suggest that at least four of these seven, 51 Pegasus, 55 Cancri, p Corona Borealis, 
and r Bootis, migrated to their present semimajor axes through the process described here. 
The planet around v Andromeda is puzzling because of its rather high eccentricity, e p = 
0.15 ± 0.04. Such an eccentricity could arise from an encounter with an object of a few M®, 
or the measured eccentricity could be in error. The last two short-period objects on the list, 
70 Virginis with M = 6.6Mj and HD 114762 with M = 10 ± lMj, are sufficiently massive 
that planetesimals scattered from their Hill sphere can be ejected with \AL\ > \AE\. If so, 
their large eccentricities could arise from the migration instability. Their relatively large a p 
might reflect the difficulty of satisfying the instability criterion (0) at small radii when M p 
is large. 

The planet orbiting 47 Ursae Majoris could have formed at its present location; however, 
we think it more likely that it migrated inward a few AU and halted because E(r) dropped 
below E c . The companion to 16 Cygni B is also likely to have migrated inward; its large e p 
could be due to interactions with the stellar companion to 16 Cygni B fllffi) . 

Equation (0) shows that a planetesimal disk mass of ~ 6 x 1O _4 M within the planetary 
orbit is required to initiate the migration instability for a planet with Jupiter's mass and 
radius, while a much larger mass of ~ 2 x 1O _2 M is required to move the planet to 0.03 
AU. Infrared observations of solar-mass T Tauri stars suggest that disks at the upper end 
of this mass range are rare, but they do occur These observations measure the total 

mass in particulate matter at r < 100 AU, with results in the range 10~ 5 — 10~ 2 M* (PZJ). 
The mass in planetesimals would be similar if the efficiency of conversion of dust particles to 



-10- 



planetesimals is high, and even higher if the observed disks already hide most of their solid 
material in planetesimal-sized objects. 

There is weak observational evidence for a correlation between the presence of short- 
period planets and high stellar metallicity We note two possible explanations for this 
correlation: (i) for a fixed mass of gas, enhancing the metallicity increases the surface density 
of planetesimals, increasing the chances for planet migration; (ii) if migration is halted by 
planetesimals hitting the star, then the planetesimals can pollute the surface layers of the 
star. 

With several planets migration becomes more complex; for example, in the solar system 
Uranus and Neptune migrate outwards rather than inwards, because Jupiter acts like an inner 
absorbing boundary similar to a nearby stellar surface (P?31). In such a system the migration 
instability is triggered at a lower surface density, since the outer planets effectively push 
the inner planet toward the star. The migration of a massive planet will reduce the surface 
density of planetesimals substantially; as a result it is unlikely that two massive planets can 
migrate to short-period orbits. It is also unlikely that any other planets in the migration 
path can survive; thus a short-period Jupiter-mass planet should have no sister planets with 
orbital radii less than a few AU. 
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Table 1. Properties of planets 



Star 


Period (days) 


a P (AU) 


Mpsmi p (Mj) e p 


t Bootis (D 


3.3128 ±0.0002 


0.0462 


3.87 


0.018 ±0.016 


51 Pegasi (TJ) 


4.229 ±0.001 


0.05 


0.47 


0.0 


v Andromedae (||) 


4.611 ±0.005 


0.057 


0.68 


0.15 ±0.04 


55 Cancri (|) 


14.648 ± 0.0009 


0.11 


0.84 


0.051 ±0.013 


p Corona Borealis (0) 


39.645 ±0.88 


0.23 


1.1 


0.028 ±0.04 


HD 114762 (D 


84.05 ±0.08 


0.3 


10 ±1 


0.25 ±0.06 


70 Virginis (@) 


116.6 


0.43 


6.6 


0.4 


16 Cygni B © 


804 


1.7 


1.5 


0.67 


47 Ursae Majoris © 


1090 


2.11 


2.39 


0.03 ±0.006 
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Fig. 1. The evolution of the i(t), a(t) and e(t) of a planetesimal having a mass 10~ 2 
that of Jupiter. The planetesimal was initially placed in the 7:4 resonance with e = 0.05 
and i = 0.05. The perturbing planet had mass \i = 10~ 3 and e p = 0.05. The planetesimal 
suffers its first close encounter with the planet after 3.5 x 10 5 yr, and is ejected after about 
3.8 x 10 5 yr. 

Fig. 2. The solid curves show the evolution of a 1 Jupiter- mass planet in a disk with 
planetesimal surface density S(r) = S Q (1 AU/r) 3 / 2 . We take initial elements as follows; the 
plane has a p = 5.2 AU and e p = 0.048, while the planetesimals have a ranging from 0.03 to 
4 AU, e in the range (e m , e m + 0.2), and i between and 0.5 radians. The labelled curves 
correspond to the following values of £© in units of g cm" 2 : A) 8000, B) 2000, C) 1200, D) 
600, and E) 40 (the nominal value for the minimum solar nebula). The corresponding gas 
surface densities in g cm~ 2 , assuming solar metallicity Z = 0.02, are A) 4 x 10 5 , B) 1 x 10 5 , 
C) 6 x 10 4 , D) 3 x 10 4 , and E) 2 x 10 3 . There is a critical value of S at ~ 200 g cm' 2 
(between D and E), below which there is very little movement. The solid circles indicate the 
positions of the various extrasolar planets (Table 1); the J indicates the location of Jupiter. 
The dotted line (labelled A*) corresponds to the case where we replaced the 2:1 resonance 
with the 3:1 resonance, in which the planetesimals are not strictly orbit-crossing unless one 
takes into account the finite Hill sphere radius. This indicates the approximate range of 
uncertainty in the final position in this idealized model. The time for the migration to occur 
is very model dependent; including more resonances will reduce the migration time. 
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